Entanglement between qubits induced by a common environment with a gap 
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We study a system of two qubits interacting with a common environment, described by a two- 
spin boson model. We demonstrate two competing roles of the environment: inducing entanglement 
between the two qubits and making them decoherent. For the environment of a single harmonic 
oscillator, if its frequency is commensurate with the induced two-qubit coupling strength, the two 
qubits could be maximally entangled and the environment could be separable. In the case of the 
environment of a bosonic bath, the gap of its spectral density function is essential to generate 
entanglement between two qubits at equilibrium and for it to be used as a quantum data bus. 
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Entangled states, showing nonlocal quantum correla- 
tions between subsystems, are indispensable in quantum 
information processing [lj. A two-qubit interaction is 
needed to implement two-qubit gates and to make en- 
tangled states. In general, there are two types of two- 
qubit couplings: direct coupling such as dipole cou- 
pling in NMR or Coulomb coupling in superconducting 
charge qubits, and indirect coupling mediated by another 
quantum system like a harmonic oscillator 0], a cavity 
mode J| > or even an electron gas 3 . A mediating system 
is expected to be separable from the qubits, so that its 
dynamics can be ignored. Recently, Cubitt et al. showed 
two qubits could be entangled by continuous interaction 
with a mediating particle in a separable state [f|. 

Entanglement of qubits and an environment causes the 
decoherence of qubits, one of the biggest obstacles in 
quantum information processing. It was believed that 
qubits should be isolated from the environment. How- 
ever, it is shown that two qubits without a direct cou- 
pling but interacting with the environment could be en- 
tangled and even in the steady state ■ Once qubits 
become entangled in the steady state, they will never be 
disentangled due to the irreversible process described in 
Ref. ■ This entangling process can not be used for the 
two-qubit coupling scheme to implement two-qubit gates. 

Like a single mediating particle, could an environment 
comprised of many particles induce the two-qubit cou- 
pling? Could it work as a quantum data bus to im- 
plement two-qubit gates at equilibrium? In this paper, 
we address this problem by studying a two-spin boson 
model shown in Fig. ^ The spin-boson model has been 
extensively studied to describe dissipative two-state sys- 
tems HQ. Here we consider two types of the environ- 
ment: a single mode case and a multi mode case. 

Single mode case. — Consider a system of two qubits 
interacting with a single harmonic oscillator with fre- 
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FIG. 1: Two circles refer to two qubits and an oval is an 
environment. The coupling between qubits and the common 
environment (solid arrow) induces the indirect interaction be- 
tween two qubits (dashed arrow). 



quency u>, described by the Hamiltonian p| 

H = tiwa) a + %X{a\ z + a^z)^ + a) , (1) 

where is a creation operator of the harmonic oscil- 
lator, ai z are the z-component of the Pauli spin matri- 
ces of the «-th qubit, and A the coupling constant. The 
Hamiltonian is solved through the canonical transforma- 
tion e x = exp[£(cri z +o- 2 z)0 + - a)} UGH- The trans- 
formed Hamiltonian reads 
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where I e g = 2hX 2 /uj is the indirect two qubit coupling 
induced by the harmonic oscillator. A total state is also 
transformed as |W) — e x \^>). 

The time evolution of the total state is given by 
= e~ x e ~ iHt / h e x |*(0)). Let us suppose an initial 
state is given by a product state |\&(0)) = \ip(0)) ® \0)e 
where \ip(0)) = <z|00) + 6|01) + c|10) + d\ll) is the ini- 
tial state of the two qubits and \0)e the ground state of 
the harmonic oscillator. Then the total state at time t is 
given by 
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where the phase 9 = I e &/h = 2X 2 /uj is due to the indi- 
rect coupling between the two qubits, the phase Tj(t) = 



2 



f— ) sinwi is due to the action of a displacement op- 
erator on the coherent state, and \a(t))E is a coherent 
state with a(t) = ^-(e~ iuJt - 1). Notice that the sub- 
space spanned by {|01),|10)} forms a decoherence free 
subspace 0. In general, Eq. © shows entanglement 
among the two qubits and the harmonic oscillator, de- 
pending on the parameters: lo and A. 

For lo ^S> A, \a(t)) is very close to |0). This means 
the overlap (0|a(t)) = exp[-r fl (i)] » 1 where T R (t) = 
(tj) 2 (1 — cosLot) goes to 0. Also the phase Ti(t) becomes 
zero. It is interesting that these two phases, Fn(t) and 
Tj(t) can also be calculated from the correlation func- 
tion 0, HG3 
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The roles of the real and imaginary parts of T(t) are 
explained below. Then, in the limit of lo 3> A Eq. JSJl can 
be written as a product state of the two qubits and the 
harmonic oscillator 



where -^(a f -a) and $(t) = e^'^^e- 
is straightforward to obtain r(i) given by 

r(t) = r B (t) + <Tj(t) 

/2A\ 2 , -| 

= 1 — I [ (1 — cos wt) — « smut \ . 



¥(*)}« o|00) + d|ll) + 
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The harmonic oscillator remains isolated always from the 
system, but it induces the indirect interaction between 
the two qubits and thus entanglement between them 5] . 
For a pure two qubits, concurrence, an entanglement 
measure, reads Cidcai = 2\ad— e l4et bc\ 0. Here the sub- 
script, ideal, stands for the case of the qubits in a pure 
state. If a = b = c= d= 1/2, we have Cidcai = | sin20t| 
as shown in Fig. 

The above analysis shows that the common harmonic 
oscillator with the high frequency, compared with the 
coupling constant to qubits, induces an inter-qubit cou- 
pling and can be used as a quantum data bus. How- 
ever, the induced qubit-qubit coupling I e g is very small, 
so that the operation time of two-qubit gates is much 
longer than that of single-qubit gates. Notice that some 
superconducting qubits use this kind of inter-qubit cou- 
pling scheme |2j. As shown below, we find that a much 
faster two-qubit gate operation is possible and the maxi- 
mally entangled state of the two qubits is obtainable even 
without the restriction of lo 3> A. 

The harmonic oscillator plays two competing roles: (i) 
making two qubits entangled by inducing the two-qubit 
coupling J c ff, (ii) having two qubits decoherent through 
entanglement with them, i.e. Tn(t). The role of the 
imaginary part T/(i) is different from that of the real part 
Tn(t). While Tn(t) makes the two qubit decoherent, i.e. 
decaying the off-diagonal elements, Tj(t) causes the two- 
qubit coupling to fluctuate as shown in Fig.|2fb). These 



could be uncovered by examining the reduced density 
matrix of the two qubits 

J26t-T*(t) 
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(7) 

which is obtained by partial tracing over the environ- 
ment, p(t) = Tr B {|*(i))(*(i)|}. Entanglement of the 
two qubits in a mixed state is measured by the concur- 
rence C = max{ri — T2 — r% — r^, 0| where r 2 are eigen- 
values of p ■ a\ y ® oiy ■ p* ■ o\ y ® a 2y • Because the total 
state is in a pure state, entanglement of the two qubits 
and the harmonic oscillator is quantified by the von Neu- 
mann entropy of the two qubits S = —Ti'(plogp). 
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FIG. 2: (color online). Concurrence C (red solid), Cidcai 
(black thick solid), entropy 2S/3 (blue dashed), and the over- 
lap exp[— Tn(t)\ (cyan dashed-dotted) as a function of time 
for lo/X = 4\/3 (a), 5 (b), and 1 (c). (d) Maximum con- 
currence C m ax, maximum entropy S majt , average concurrence 
Cavg j and average entropy SVvg as a function of n = (a>/4A) 
for < 6t < it/2. 

As lo/X becomes large, C goes to Cidcai and S to 
as depicted in Fig. |3 This is explained by means of 
the Born-Oppenheimer approximation that is based on 
the assumption of the weak coupling between the two 
qubits and the environment. The frequency lo of the har- 
monic oscillator is larger than that of the two qubits, 9. 
The harmonic oscillator oscillates very fast in comparison 
with the two qubits. So the two qubits feel the harmonic 
oscillator stays in the same state. However, the condition 
of lo 3> A is not a unique way to make the two qubits en- 
tangled maximally. Surprisingly, we find the maximum 
entanglement of the two qubits, C = 1, at 6t = tt/4 un- 
der the condition that lo/X = 4y / n with n = 1,2, . . . as 
shown in Figs.Efa) and^d). This is due to the fact that 
the frequency 9 for entangling the two qubits is commen- 
surate with the frequency lo of the harmonic oscillator. 
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If the condition of cj/A = 4y / n, n = 1, 2, . . . is not met, 
then the concurrence C does not reach 1 at 9t = ir/4 
and at 9t = n/2 due to T(t) as shown in Fig. |2fb). 
The oscillation period of C (red solid line) does not coin- 
cide with that of the ideal case Cideai (thick black line). 
This implies that the two qubit coupling fluctuates due 
to Ti(t). For lo/X < 1, our numerical study shows that 
the two qubits could not be entangled. Fig. (2fc) shows 
the case of oj = A. 

Fig. Efd) shows two competing roles of the harmonic 
oscillator. Let us define the average concurrence C avg 
for a period r = n/29 by C avg = \ J T C(t)dt. For the 
first period < 9t < tt/2, the maximum concurrence 
Cmax reaches the peak at positive integers n and the av- 
erage concurrence C avg increases as w/A increases. In 
contrast, the maximum entropy S niax and the average 
entropy S avg = - J S(t) dt decrease. Due to the corn- 
mensuration of two frequencies u and 9, S avg and C avg 
show the behavior of the stair case. 

Two-spin boson model. — In the single mode case, the 
dynamics is reversible. If the single harmonic oscillator 
is replaced by the bosonic bath, the dynamics of qubits 
becomes irreversible in the sense that a pure initial qubit 
state becomes mixed and never returns to a pure state. 
The total Hamiltonian of the two-spin boson model is 
given by 



N 



N 



(8) 



where the first term describes the bosonic bath and the 
second the interaction between the two qubits and the 
bosonic bath (environment). Here it is assumed that N 
goes to infinity. The Hamiltonian of two qubits H q — 
J2i=i.2 (BixVix + B iz a iz ) is set to be zero by turning off 
an external magnetic field. Also it is assumed that no 
direct two qubit coupling is available. 

The environment is characterized by the spectral den- 
sity function J(uo) = . \? S(uj — u>j) As shown in 
the inset of Fig.|3J-(a), let us consider an Ohmic environ- 
ment with a gap ojq and an exponential cutoff function 
of the cutoff frequency w c 



J{lo) = a(u> - w ) e -(^-^o)M 0( w . 



wo) , (9) 
a step func- 



where a is the coupling parameter and 9(x) 
tion. 

Along the same lines of the single mode case, the 
canonical transformation 



exp 



{<7z 



3 J 



(10) 



eliminates the coupling between the two qubits and the 
common environment in the Hamiltonian (JSJ and induces 
the indirect coupling between the two qubits. The trans- 
formed Hamiltonian becomes 



H = ^ ftiOjCtjClj - icffCi z O- 2 z - IcS ■ 



(11) 
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FIG. 3: (color online). (a) exp[— Tr^)] as a function 
of cuct and (b) its log scale plot. Legends f),...,6) re- 
fer to (u> /u c ,a) = (0,0.25), (0.01,0.25), (0.1,0.25), (0,0.5), 
(0.01,0.5), (0.1,0.5), respectively. Inset of (a): the spectral 
density function J(uj) with a gap ojq. (c) Concurrence C 
(solid), entropy 2S/3 (dotted), and exp[— Ya(t)] (dashed) as 
a function of Qt. 



where the induced coupling is given by 



Jeff = %T— = 2h 

* * 1,1; 



w 



diu . 



(12) 



For the Ohmic environment without a gap, ujq — 0, we 
have I ff = 2hauj c . 

Let us investigate the dynamics of entanglement. For 
an initial state |*(0)) = |^(0)) (g) \0} E where \0) E = 
|0i, O2, ■ ■ ■ , On) is the ground state of N harmonic oscil- 
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lators, the state at time t is given by 



!*(*)) 



o|00) ® |a(t)) B + rf|ll> <8> |-a(t))j 



e- aet+lT ' {t \b\Ql) + c|10))<8 |0) 



(13) 



where \a(t)) E = |{^(e" ia,it - 1)}) is the collection of 
the coherent states. The imaginary part of T(t) becomes 



r /(t )=£'' 2A ' 



sin Wjt 



— — sin (art) aw . 



(14a) 
(14b) 



Also we obtain the overlap B(0|a(t))B = exp[— r#(i)] 
where 



r*(t) = 4 / 

Jo 



£M (1 



cos c*;t) dcj . 



(15) 



From Eq. , the density matrix of the two qubits be- 
comes the same form of Eq. (JTJl . For the Ohmic spectral 
density @ without a gap, uiq = 0, we have 0, 0] 



r fl (f) = 2aln[l + (cj c tf 
r 7 (t) = 4atan~ 1 (cj c £) . 



(16a) 
(16b) 



The slopes of solid lines in Fig. 01b) are 4a. In this 
case exp[— Tn(t)] goes to zero as time goes on, but never 
returns to 1 in contrast to the single-mode case. Then 
the density matrix of the two qubits eventually becomes 
a mixed state with the maximum entropy S — 2/3. For 
the environment without a gap, the two qubits become 
entangled and disentangled only for the transient times 
not in the "steady state" as shown in Fig.^c). Here the 
term "steady state" is used in the sense that in the steady 
state the entropy S(t) of the two qubits does not change 
in time 0] . In this case, there is no further information 
flow from the two qubit to the environment. 

It is clear that in the steady state the two qubits could 
be entangled and disentangled only if exp[— r#(t)] re- 
mains finite. How is it possible? We find that the 
small gap loq <C u> c of the Ohmic environment makes 
exp[— Tn(t)] finite in the steady state and the two qubits 
entangled. Fig. [3] shows that for the Ohmic environment 
with a gap the decohcrcnce term exp[— r#(f)] does not 
vanish in the steady state. The role of the gap is to pre- 
vent the complete drain of information from the qubits 
to the bath. We would like to mention a similar phe- 
nomenon, that is, the inhibition of the spontaneous emis- 
sion of an atom in photonic-band-gap materials . No- 
tice that the gap of the spectral density function causes 
the phenomenon called false decoherence 16]. It would 
be interesting to derive a master equation of the two 
qubits in the Lindblad form when the spectral density 
of the environment has a gap. 




FIG. 4: (color online), (a) Maximum concurrence of the two 
qubits (b) entropy of two qubits at equilibrium as a function of 
a and u>q/u> c . (c) exp[— r_R(txj)] as a function of temperature 
T and u>o/^c for a = 0.25. 



Figs. Ufa) and Efb) show the maximum concurrence 
C max and the entropy S in the steady state at zero tem- 
perature. C max (S) decreases (increases) as a increases 
even for ujq ^ 0. The effect of the non zero temperature 
of the environment is shown in Fig.^Jc). For finite tem- 
peratures, we take the initial density matrix W(0) of the 
total system given by 

w(o) = \m)(m\®PE, pe = e HE J kBT , (17) 

where the environment is in a thermal equilibrium state 
Pe at temperature T. Here He — fiu>ia\a,i and Ze 
is the partition function of the environment. At finite 
temperatures, the decoherence factor is given by |3.l3.ITo| 



r«(t) 



J(w) 



cotMsffrX 1 -cosLut)(ku (18) 



which decreases the effect of the gap as shown in Fig.0Jc). 
However, the real part Tn(t) is identical to Eq. ((HI) . 
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We would like to make some remarks. First, it is easy 
to obtain the induced interaction for many qubits inter- 
acting with a common environment. It has the form of 
pair- wise interactions H oc I c s < J izO~j Z - In this paper 
no external magnetic fields are applied on qubits, i.e., 
Bi, i — 1,2. For 7^ 0, one can also obtain the effec- 
tive Hamiltonian using the flow equation renormalization 
group 01 1 which is a generalized canonical transforma- 
tion. 

Second, there is an interesting physical systems similar 
to the two-spin boson model. Ref. [5[ considers two spins 
of the two quantum dots coupled with a 2-dimensional 
electron gas. In this case the induced interaction is noth- 
ing but the Ruderman-Kittel-Kasuya-Yosida (RKKY) in- 
teraction. However, the spin on a quantum dot becomes 
entangled with a 2-dimensional electron gas 0]. Also 
the density of states of an electron gas has no gap. Ac- 
cording to our theory in this paper, although the two 
spins are coupled via the RKKY interaction, their state 
is not pure but mixed. A further study should be done 
in order to know whether the RKKY interaction could 
be used for a two-qubit coupling scheme in solid state 
quantum computers |^. 

In conclusion, we have studied the two qubits interact- 
ing with a common environment with a gap, described 



by a two-spin-boson model. For the environment of a 
single harmonic oscillator, if the induced two-qubit cou- 
pling strength is commensurate with the frequency of the 
harmonic oscillator, the two qubits could become maxi- 
mally entangled. Also we discussed the roles of the real 
and imaginary parts of T(t). In the case of the two-spin 
boson model with the spectral density function of the 
gap ujq, at the steady state the two qubits become en- 
tangled and disentangled due to the indirect interaction 
between qubits induced by the environment. Our result 
shows explicitly that two qubits can be entangled due to 
the environment with a gap while the environment itself 
is not entangled with qubits. 
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